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>■ ■ A joint spectral analysis of some Chandra ACTS X-ray data and Molonglo Ob- 



servatory Synthesis Telescope radio data was performed for 13 small regions along 
C^^ ' the bright northeastern rim of the supernova remnant SN 1006. These data were 

fitted with a synchrotron radiation model. The nonthermal electron spectrum 
(^ ■ used to compute the photon emission spectra is the traditional exponentially cut 



off power law, with one notable difference: The power-law index is not a constant. 
It is a linear function of the logarithm of the momentum. This functional form 
>. ■ enables us to show, for the first time, that the synchrotron spectrum of SN 1006 

j_| ' seems to flatten with increasing energy. The effective power-law index of the elec- 

tron spectrum is 2.2 at 1 GeV (i.e., radio synchrotron-emitting momenta) and 
2.0 at about 10 TeV (i.e.. X-ray synchrotron-emitting momenta). This amount 
of change in the index is qualitatively consistent with theoretical models of the 
amount of curvature in the proton spectrum of the remnant. The evidence of 
spectral curvature implies that cosmic rays are dynamically important instead 
of being "test" particles. The spectral analysis also provides a means of deter- 
mining the critical frequency of the synchrotron spectrum associated with the 
highest-energy electrons. The critical frequency seems to vary along the north- 
eastern rim, with a maximum value of l.l^o 5 ^ ^O^^ Hz. This value implies that 
the electron diffusion coefficient can be no larger than a factor of ~4.5-21 times 
the Bohm diffusion coefficient if the velocity of the forward shock is in the range 
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2300-5000 km s~^. Since the coefficient is close to the Bohm hmit, electrons are 
accelerated nearly as fast as possible in the regions where the critical frequency 
is about 10^^ Hz. 

Subject headings: acceleration of particles — cosmic rays — ISM: individual 
(SN 1006) — radiation mechanisms: nonthermal — supernova remnants — X- 
rays: general 



Introduction 



Several observational and theoretica l clues support the sugge stion that Galactic cos- 
mic rays, up to an energy of 100 TeV (JLaeage fc CesarskyI Il983l ) or more ( Jokipiil 119871 : 
Volk fc BiermamjIlQSSi : iBell &: Lucekll200ll ). are accelerated predominantly in the shocks of 
supernova remnants. If Galactic supernovae occur at an average rate of one event every 30 
years, then an average supernova remnant must transfer about 10% (JDrury et al.lll989l ) of 
the initial 10^^ ergs of kinetic energy of the ejecta to cosmic rays. In this case, the cosmic-ray 
energy density may be large enough to affect the structure of the shock (i.e., cosmic rays 
may not be mere "test" particles). 

Three consequences of a large cosmic-ray pressure are potentially observable. One con- 
sequence is that the ambient material is s lowed before it crosses the subshock. This effect re- 
duces the temperature of the shocked gas (IChevalierlll983l : Ellison et al.ll2000l : lDecourchelle et al 
20001 ). An upper limit on the temperature is provided by the well-known relation between 
the shock speed Vg and the postshock temperature of a gas whose ratio of specific heats 
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kTi < — miV^ 
16 ' 



(1) 



where k is Boltzmann's constant and rrii and Tj are the mass and immediate postshock 
temperature, respectively, of particle species i. Since kinetic energy is transferred to cosmic 
rays at the expense of the thermalization of the shocked gas, t he equality is approp riate 
only in the limit that the cosmic-ray energy density is negligible. iHughes et al.l (J2000h ana- 
lyzed the transverse motion of the X-ray-emitting material (mostly reverse-shocked oxygen 
and neon) in the supernova remnant IE 0102.2—7219 and inferred a forward-shock velocity 
Vs = 6200^1600 km s~^. In this case, the temperature of the shocked protons must be lower 
than the right-hand side of equation ([1]) or else the fitted el ectron temperature will be too 
low to be explained by Coulomb heating (JHughes et al.l 120001 ). However, the proton temper- 
ature is sensitive to the shock velocity (T oc f^), and a variety of speeds are reported for 
IE 0102.2—7219. iFlanagan et al.l (120041 ) report that the Chandra High Energy Transmission 
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Grating Ne X line emission data are consistent with a na odel that includes radial Ne X ve- 
locities up to 1800 ± 450 km s~^. iFinkelstein et al.l (120061 ) analyzed several fe atures in O III 



i mage s. The mean transverse velocity of the features is 2000 ± 200 km s~^. lEriksen et al. 
(I2OOII ) report that the motion of some O Ill-emitting material is best described by a radial 
O III velocity of about 1800 km s~^. Since these velocity results differ and since none of 
the measurements provides a direct measure of the velocity of the forward shock, the claim 
that a significant fraction of the internal energy in IE 0102.2—7219 has been transferred to 
cosmic rays requires additional support. 



A second consequence is tha t the total compression ratio is larger than 4 (JEUison &: Reynolds 



I991I ; iBerezhko fc EUisonl Il999l ). Since cosmic rays slow the upstream material before it 
crosses the subshock, the velocity of the shocked material relative to the subshock is re- 
duced. As a result, the contact discontinuity is closer to t h e sub shock than it would be 
in the absence of a large cosmic-ray pressure. I Warren et al.l (120051 ) report that in Tycho's 
supernova remnant the mean ratio of the radius of the contact discontinuity to the radius of 
the forward shock is .93=1=0.02. This separation corresponds to a total compression ratio 



1+1.9 



r = 5.1]'^i;o. Similarly, ICassam-Chenai et al.l (l2008l ) report that for the southeastern rim of 



SN 1006 the ratio varies from a value of 0.911qq2 (i-^-; '^ 



4.0^1'^) between the bright X-ray 



synchrotron-emitting filaments to a value of 1 (i.e., r = 00) along the filaments. The mean 
value between the filaments is 0.96 ± 0.03 (i.e., r = 9.0]'^3^g). These authors explore ways 
in which the location of the contact discontinuity can approach the location of the forward 
shock, but they cannot explain why these two features curiou sly appear to be coincident 
with one another along the filaments. iKsenofontov et al.l (120051 ) . who also studied SN 1006, 
infer a mean compression ratio r = 5.2. 

A third consequence of a la rge cosmic-ray pressure is that the shock transition region 
is broadened or "smeared out" (JEllison &: Reynolds! Il99ll : iBerezhko fc EUisonl Il999l ). Only 
the subshock has a short transition length. In this case, low-energy cosmic rays, which have 
relatively small diffusion lengths, experience only a portion of the velocity gradient as they 
scatter back and forth across the subshock. Higher energy particles, which have relatively 
large diffusion lengths, experience a larger portion (or all) of the velocity jump. Since the 
rate of energy gain increases as the velocity difference increases, higher energy particles gain 
energy faster than lower energy particles. As a result, cosmic- ray spectra do not have power- 
law distributions. The spectra fl atten with in c reasin g ene rev (lBelllll987l:lEllison fc Reynolds 
I991I : IBerezhko fc EUisonl 1 1999h . Ijones et all J2003h and Ivinkelalj (|2006|) report evidence 
of curvature in the synchrotron spectra of Gas A and RGW 86, respectively. One possible 
explanation for such curvature is that the cosmic-ray electrons producing the synchrotron 
emission have curved spectra. 



-4- 



Here we describe the results of an analysis of some radio and X-ray synchrotron data, 
which suggest that the synchrotron spectrum of SN 1006 might be curved. The data, as- 
sumptions, and analysis techniques are described in § [21 The results of the analysis are 
discussed in § [31 Our conclusions are summarized in § [H 



2. Data and analysis 

X-ray spectra of the filaments along the northeastern rim of SN 1006 were obtained by 
analyzing 68 ks of data from the 2000 July 10-11 Chandrsu observation of the source. The 



result s of p revious analys e s of t he data from this observation were published by [Long et al. 



(I2OO3I ) and iBamba et al.l (120031 ). During the observation, the telescope was pointed at a 
location near the middle of the northeastern rim (a = 15^03™51.56'^, 6 = — 41°51'18.8"; 
J2000). The X-rays were detected using the Advanced CCD Imaging Spectrometer (ACIS), 
which includes a 2 x 2 (ACIS-I) and a 1 x 6 (ACIS-S) array of CCDs. Six of these 10 detectors 
(ACIS-I2, -13, -SI, -S2, -S3, and -S4) were used for the observation. Each 1024 pixel x 
1024 pixel CCD has a field of view of 8.4' x 8.4'. The angular resolution of the Chandra 
mirrors and ACIS varies over the observed portion of SN 1006 from about 0.5" at the aim 
point to 20" for a region that is 17.5' off-axis. The on-axis effective area for the mirrors and 
ACIS-S3 has a maximum of about 720 cm^ at 1.5 keV and is greater than 10% of this value 
for energies between about 0.3 and 7.3 keV. The fractional energy resolution {FWH.M/ E) 
between these energies ranges from about 0.4 to 0.03, respectively. The sensitive energy 
bands and energy resolutions of the other five CCDs used for SN 1006 are typically narrower 
and worse, respectively, than the energy band and resolution of ACIS-S3. 

The ACIS data were filtered to remove the events that (1) have GRADE = 1, 5, or 7, 
(2) have one or more of the STATUS bits set to 1 (except for events that have only one or 
more of the four cosmic-ray "afterglow" bits set), (3) occur on a bad pixel or column, (4) are 
part of a serial readout streak on ACIS-S4, or (5) occur in the time interval during which 
the mean background count rate was more than twice the nominal rate (i.e., frames 20515- 
21294). An image of the 66 ks of "good" data is displayed in Figure [H The 13, 100 pixel x 
100 pixel (49" x 49") boxes are the regions of the bright, northeastern filaments used for the 
X-ray spectral analysis. The PHA spectra, ARFs, and RMFs for each region were created 
using version 2.3 of the CIAOg tools dmextract, mkarf , and mkrmf , respectively. 



^For more information about the Chandra X-ray Observatory refer to the Chandra Proposers' Observatory 
Guide, which is available at http://asc.harvard.edu/proposer/POG/index. html, 

^For more information about the CIAO analysis tools, refer to http://asc.harvard.edu/ciao 
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The X-ray analysis was performed using data in the energy band 2-7 keV. The higher 
energy data were excluded because they are dominated by background events. The lower 
energy data were excluded for three reasons. First, a small, but noticeable, amount of thermal 
emission is evident at the lower energies. Therefore, inclusion of the lower energy data 
requires a spectral model that includes a thermal emission component. The extra parameters 
required to fit the thermal emission make the spectral fits unnecessarily co mplicated. By 



ignor ing the data below 2 keV, nearly all of the thermal emission is excluded (JBamba et al. 



20031 ). As displayed in Figured] (and as originally reported by Koyama et al. 1995), the 
high-energy nonthermal emission is concentrated along the rim. The second reason is that 
the ACIS detectors are not well calibrated at energies below about 0.7 keV. The last reason 
is that the absorption column density can be frozen at some preferred value to simplify the 
fitting pr ocess. Our pref erred value is ri} 



6.0 X 10^*^ cm ^, which is consistent with the 



results of lSchaeferl fll996[ ) and lAUen et all fl200l[ ). 



The radio data used for the spectral analysis include the compilation of flux density 
measurements listed in Table [H These measurements are based on analyses of the emission 
from the entire supernova remnant. In the cases where no flux density uncertainties are 
reported, the uncertainties are assumed to be 10% of the reported flux densities. The results 
of the spectral analysis are insensitive to the actual fraction used at least for fractions in 
the range 5%-20%. The values hsted in Table [1] are plotted in Figure [2j The dashed line 
in this figure is the result obtained when only the radio data are fitted with a power-law 
model: S{iy) = 17.6t^j[z//(l GHz)]"" Jy, where a = 0.60+J]:|]^. Note that the uncertainties 
for the power-law index and normalization, which are quoted at the 90% confidence level, 
are not indepe n dent. This result is an update of and consistent with the result reported 
bv lAUen et aP J200lh : S{u) = 17.9 ± 1.1^(1 GHz)]"" Jy, where a = 0.57 ± 0.06. The 



dotted line in Figure [2] is a power-law model with the parameters reported by iGreenI (120061 ): 

S{u) = 19[z//(l GHz)]-o-6 Jy. 

The radio spectrum for the entire remnant is used to construct radio spectra for each 
region because th ere is little evidenc e of ra dio sp ectral variability. For exa mple, a comparison 
of the images of IStephenson et al.l (119771 ) and [Gardner fc Milnd (119651 ) suggests that the 
slopes of the 0.408-2.7 GHz radio spectra of the northeastern and southwestern rims are 
about the same. It is possible that the slope may vary along each rim, in which case 
the best-fit amounts of curvature for individual regions may be inaccurate, but the mean 
amount of curvature for the 13 regions is relatively insensitive to such variations. Two sets 
of radio data were used with each X-ray spectrum. Both sets are obtained by multiplying 
the remnant-integrated flux densities and uncertainties of Tabled] by scaling factors (. One 
set of factors (Tables [2] and [3|) is obtained by dividing the 843 MHz flux densities for the 
13 regions marked with black squares in Figure [3l from the Molonglo Observatory Synthesis 
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Telescope (MOST), by the 843 MHz flux density for tlie entire remnant. Tliese regions are 
tlie same ones used to obtain tlie X-ray spectra (Fig. [1]). Tfie second set of factors (Tables H] 
and E]) is obtained in the same way, except that the flux densities for the 13 red squares 
are used instead of the flux densities for the black squares. The red squares are the regions 
where the flux densities peak. That is, the flux density along a line from the center of 
the remnant through the center of a black square peaks at the location of the center of 
the corresponding red square. Although both sets of flux densities represent local averages 
because the half-power beamwidth associated with the 843 MHz MOST image (44" x 66"; 
see Fig. [3]) is comparable to the size of the regions (49" x 49"), they are reasonable estimates 
of the minimum (black) and maximum (red) flux densities for the 13 regions. Therefore, 
using these two sets of radio spectra provides a means of obtaining upper (Table [2]) and 
lower (Table H]) limits on the amount of curvature for each region. 

Version 1.4.8 of the spectral fitting package ISIqfl (JHouck fc DeNicolall200Cll ) was used to 



fit the data for each region with a model that includes a synchrotron radiation component 
for the X-ray and radio data and an interstellar absorption component for the X-ray data 
(the XSPEC model wabs, with the relative abundances of Anders & Grevesse 1989). Since 
the 2-7 keV emission is dominated by nonthermal emission, no thermal X-ray component is 
included. 

The synchrotron spectra are based on a nonthermal electron spectrum of the form 



^ = A f ^V'^"'"^"^xp f^^^ L (2) 

dp \PoJ \ Pm J 

where n is the electron number density, p = 'jmv (7, m, and v are the Lorentz factor, rest 
mass, and velocity of a particle, respectively), A is the number density at p = po (in units 
of Po ^ cm"'^), po = 1 GeV c~^ (c is the speed of light), F is the differential spectral index 
at p = po, a is the spectral "curvature," the logarithm is base 10, and Pm is the exponential 
cutoff (or "maximum") momentum. This form is the same as the standard power law with 
an exponential cutoff, except for the term a log (p/po)- Introduction of this term produces an 
effective spectral index, Fcs = F — a log (p/po), that is a linear function of the logarithm of 
the momentum. For example, at momenta p of lO'^po; lO^po^ lO^po; and IO^Pq, the effective 
differential spectral indices are F, F — a, F — 2a, and F — 3a, respectively. If a > 0, then 
the spectrum flattens with increasing momentum. If a < 0, then the spectrum steepens 
with increasing momentum. If a = 0, then the spectrum has no curvature and equation ([2]) 
reduces to the standard power law with an exponential cutoff. An exponential cutoff is 



'For more information about ISIS, see http://space.mit.edu/cxc/isis 
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used in equation (121) fo r sim plicity. Note that lEUison et al.l (120011) , lUchiyama et al.l (l2003l ) 



Lazendic et al.l (120041 ). and 



Zirakashvili fc Aharonianl (120071 ) use a mo r e gen eral form for 



the cutoff (exp[— (p/pm)"], where s = |, ^, 1, or 2) and that iProtherod (120041 ) performed a 
detailed study of the shape of the cutoff for particles experiencing diffusive shock acceleration. 

The synchrotron spectrum is computed by evaluating the following integral expression 
to obtain the differential photon flujQ (in SI units): 



dF Vs V^e^B 



d(hiy) And'^ 2eQmch'^iy 



dp 



dn 
dp 



de sin^^ 



V 



U/Vc 



dxK5/3 (x) 



(3) 



where Vs is the volume of the synchrotron-emitting region, d is the distance of the source 
from Earth, e is the unit of electric charge, B is the total magnetic field strength (not B±), 
eo is the permittivity of free space, h is Planck's constant, z/ is the frequency of an emitted 
photon, 6 is the pitch angle between the electron momentum and magnetic field vectors, 
7^5/3 is an irregular modified Bessel function of the second kind, and 

3e 



{p^c^ + m^c'^) B sine 



(4) 



Atttti^c^ 

is the critical frequency of the synchrotron spectrum produced by an electron with a mo- 
mentum p. The frequency at which a synchrotron power spectrum peaks is Up = 0.286z/c- 
Hereafter, only the critical frequency Uc is used. 

We refer to the critical frequency associated with electrons that have a momentum equal 
to Pm as the "cutoff critical frequency" or just the cutoff frequency to emphasize that it is 
associated with electrons at the cutoff momentum of the electron spectrum [z/^ = i^c{p = 
Pm)]- The distribution of pitch angles 9 is expected to be nearly isotropic because the 
forward-shock speed (2300-5000 km s~^; Laming et al. 1996; Dwarkadas & Chevalier 1998; 
Ghavamian et al. 2002) is much less than the speed of light. In this case, the mean value of 
sin 6 in equation (]!]) is about 7r/4, 



Urn = 1.26 X 10 
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Pn 



10 TeV c- 



2 / B 

iooTIgJ ^^' 



and 



Pr^ 



^"- = °-'''UOTeVc-; VlOO^G 



B 



keV. 



(5) 



(6) 



During the spectral fitting process, the right-hand side of equation (^ is computed 
each time the spectral fitting program changes the parameters of the model. The complete 



*For more information about this model, see lHouck fc Allen! ( 20061) . 



set of spectral parameters includes the spectral index F, the curvature parameter a, the 
cutoff momentum p^, the total magnetic field strength B, the normalization constant A^^ = 
AVs/'ind'^, and the absorption column density rin. The spectral index is determined by the 
radio and X-ray data. The curvature is sensitive to the relative radio and X-ray fluxes. Since 
only two of the three parameters Pm, B, and Ns are independent, pm was fixed and the cutoff 
frequency {um oc p^B] eq. [5]) is reported here instead of the magnetic field strength. This 
frequency is sensitive to the shape of the X-ray spectrum. The same normalization is used 
for both the X-ray and radio data for a given region because the nonthermal X-ray and radio 
emission from the region is assumed to be produced by a common population of electrons. 
The normalization A^^ is determined b y the radio data. The absorption column density rin 



was fixed at the value 6.0 x 10^° cm'^ flSchaefedll996l : lAUen et al.ll2001f ) 



Four spectral fits were performed for each region. The first pair are identical to one 
another except that the curvature parameter is a free parameter for one (Table [2]) and is 
fixed at a value of zero (i.e., no spectral curvature) for the other (Table [3]). In both cases, 
the radio spectra are normalized to the flux densities in the regions that are cospatial with 
the X-ray regions (i.e., the regions marked with black squares in Fig. [3]). Since the cospatial 
radio spectra represent lower limits on the radio fluxes, the amounts of curvature listed in 
Table [2] are upper limits. The second pair of fits (Tables H] and [5]) is identical to the first 
pair except that the radio spectra are normalized to the flux densities in the "peak" regions 
(i.e., the regions marked with red squares in Fig. [3l) instead of the cospatial regions. Since 
the peak radio spectra represent upper limits on the radio fluxes, the amounts of curvature 
listed in Table H] are lower limits. For each of the four sets of fits, the best-fit values of the 
spectral index vary little from region to region, because the same radio spectral shape was 
used for each region. The best-fit model for region 6, the region with the largest number of 
counts, is plotted in Figures [2] and HI 



3. Discussion 

It is not immediately obvious whether the cospatial or peak radio fluxes yield results 
that are more accurate, because the X-ray and radio synchrotron morphologies are not the 
same. For example, along 49"-wide strips that pass through the 13 X-ray regions toward the 
center of the remnant, one finds that the X-ray emission peaks at the locations of the black 
boxes in Figure [3] and that the radio emission peaks at the locations of the red boxes. There 
are at least three possible explanations for this difference. One explanation is that the TeV 
electrons that produce the X-ray synchrotron emission suffer significant synchrotron losses 
as they diffuse downstream from the shock. These losses deplete the highest-energy end of 
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the electron spectrum, causing the X-ray synchrotron emission to dechne in the downstream 
region. As a result, the X-ray peak is relatively narrow and confined to a region near the 
shock. Since synchrotron losses are negligible for the GeV electrons that produce the radio 
emission, the radio emission is not attenuated in the same manner. The radio peaks can be 
broader and farther from the shock than the X-ray peaks. Nevertheless, there is no reason, 
on these grounds, to use different radio and X-ray extraction regions, because the radio and 
X-ray synchrotron emission for a region is produced by the same population of electrons. 
The radio emission is produced by the GeV portion of the electron spectrum and the X-ray 
emission is produced by the TeV segment. Synchrotron losses do not introduce a spatial 
separation between the GeV and TeV electrons. 

Another explanation for the different X-ray and radio morphologies does involve an 
apparent (not actual) spatial separation. The X-ray and radio instruments do not yield 
images that have the same spatial resolution. The size of the spatially dependent point-spread 
function of the Chandra X-ray telescope is much smaller than the size of the 49" x 49" boxes 
and, hence, is neglected. However, the size of the beamwidth for the MOST image (44" x 66") 
is not negligible, because it is comparable to the size of the boxes. The implications of the 
difference in the qualities of the X-ray and radio images are illustrated in Figure O The black 
histogram in the top panel shows the radial X-ray profile for region 6. The red curve in this 
panel is obtained by using a Gaussian function to smooth the histogram to the beamwidth 
of the radio instrument. Note that the peak of the red curve is about 18" downstream from 
the peak of the histogram, which suggests that the radio extraction region should be shifted 
toward the radio peak by the same amount. Similarly small shifts are obtained for the other 
12 regions and for a simple model of a uniformly emitting spherical shell (Fig. [5l bottom). 

A third explanation is that the bright radio-emitting region has expanded in the 17 
intervening years between the time when the MOST data were obtained and the time of 
the Chandra observ ation. Using the mean radio expansion result of 0.049% ± 0.014% yr~^ 



(iMoffett et al.lll993l ) suggests that the radio peak has shifted about 7.6" ± 2.2". To com- 
pensate for this expansion, the radio extraction region should be shifted toward the radio 
peak by the same amount. Since this shift and the apparent shift described above are small 
compared with the size of the extraction region, the following discussion focuses primarily 
on the results obtained using the cospatial radio fluxes. 



3.1. Curvature 

While the values of x^ per degree of freedom (dof ) are acceptable for each of the four 
sets of fits, a comparison of Tables [2] and [3] (and of Tables H] and E]) reveals that the values 
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are uniformly lower when curvature is included as a free parameter. For example, Figure [6] 
shows the 1, 2, and 3 a confidence contours for region 6 in the parameter space defined by 
the electron spectral index F and the spectral curvature a. If the cospatial radio flux density 
is used, the results for this region suggest that an uncurved spectrum (i.e., the dashed line 
at a = 0) can be excluded at about the 2.7 a confidence level. Similar results, with varying 
degrees of statistical significance, are obtained for the other 12 regions. Table [2] lists the 
probability -Pax^ that the reduction in the value of x^ when curvature is included as a free 
parameter is due to chance. The majority of the probabilities are less than 1%. Some are 
much less. The probabilities obtained using the F-test (Pf) are within a factor of 2 of the 
values listed in Table [2l except for regions 1, 2, 3, 5, 6, and 7, where Pp = 0.058, 0.0035, 
9.9 X 10~^, 0.0013, 0.0030, and 0.048, respectively. The relatively large discrepancies between 
Pax^ and Pp for these regions arise because the F-test depends on the value of x^/dof in 
addition to the value of Ax^. Note that the values of x^/dof are significantly less than 1 
for regions 1, 2, and 3 and that regions 5, 6, and 7 have the largest values of x^/dof. Since 
the values of Pax^ are typically larger (i.e., less significant) than the values of Pp, only the 
values of Pax^ are listed in Table [2l Collectively, these probabilities show that the spectra 
are significantly better fitted by a curved synchrotron model than an uncurved one. 

The strength of this evidence depends on the location used to obtain the radio flux 
density. The radio flux density for a region increases from the cospatial location to the peak 
location. Since the appropriate location to use lies between these two, the cospatial and peak 
flux densities represent the lower and upper limits, respectively, on the radio flux density for 
the region. The amount of curvature depends on the relative X-ray and radio synchrotron 
fluxes. Therefore, the relatively small cospatial flux density yields an upper limit on the 
amount of curvature and the peak flux density yields a lower limit. Although the curvature 
values in Table H] are lower limits, even these results favor a curved synchrotron model. Yet, 
as noted above, the upper limit values in Table [2] are expected to be closer to the actual 
amount of curvature because the appropriate location to use for the radio flux density is 
much closer to the cospatial location than the peak location. 

The evidence of curvature cannot be attributed to an instrumental effect. It must be 
associated with the synchrotron spectrum of the source. For example, the need for curvature 
could be ehminated if the effective area of the X-ray detector system were too low, but the 
area would have to be in error by an implausibly large factor of about 4.6. The use of a 
curved source model is also compelling in the sense that the values of the best-fit indices 
and curvature parameters are consistent with the expected values. One expectation is that 
the best-fit index should be consistent with the radio spectral index. That is, the best-fit 
model should fit the radio spectrum instead of simply intersecting it. The mean value of 
the index F = 2.221 ± 0.013 if curvature is included as a free parameter (Table [2]). This 
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value is consistent with the index obtained if only the radio data are fitted with a power- 
law model (r = 2q; + 1 = 2.20lo;Jg; -^^S- 12)- Yet this result is not surprising, because the 
fitted value of the index is determined largely by the shape of the radio spectrum. For 
comparison, the mean value of the index F = 2.031 ± 0.017 if an uncurved model is used 
(Table [3]). This value differs significantly from 2.22 because the uncurved index, unlike the 
curved index, is sensitive to the relative radio and X-ray fluxes. While this result seems to be 
closer to the case of having the model simply intersect the radio spectrum, an index of 2.03 
is still marginally consistent with the radio data because the radio data have relatively large 
uncertainties. Otherwise, the radio data might provide a means of discriminating between a 
curved and an uncurved model. 

Another expectation is that the value of the curvature parameter should be positive 



That is, the spectrum should fl atten with increasing energy (JBelll 119871 : lEUison fc Reynolds 



199ll : iBerezhko fc Ellisorull999l ). The results in Table [2] show that the best-fit values are 
greater than zero for all 13 regions. If the spectra are not curved, then the use of a curvature 
parameter may not be physically meaningful and one might expect roughly half of the values 
to be negative. If either sign is equally likely, then the chance probability that all 13 values 
have the same sign is 0.00024 (2~^^). 

A third expectation is that the highest-energy electrons, the ones that experience the 
full strength of the shock an d that produce the X- ray synchrotron emission, should have an 



effective spectral index of 2 (IBerezhko et al.ll2002l ). If equation ([2]) is an accurate represen- 



tation of the shape of the electron spectrum, then the effective spectral index in the X-ray 
band Tx = F — alog (px/po), where px is the fiux-weighted mean momentum of the electrons 
that produce the X-ray synchrotron radiation. To compute px from the fitted value of Um 
requires knowledge of the magnetic field strength (eq. [5]). Since the results of our analysis 
do not yield an estimate of the field strength, we assume that px = 10 TeV c~^. (The results 
are rather insensitive to this assumption as long as px is within an order of magnitude of 
10 TeV c~^.) In this case, Tx = 2.005 ± 0.027. This index and the index obtained using the 
uncurved model (F = 2.03l;!^o;Qjg) are both consistent with the idea that the highest-energy 
electrons in SN 1006 have a spectral index of 2.0. 

If the effective spectral index of the X-ray synchrotron-producing electrons Tx = 2 
(aside from the exponential cutoff), then the expression for the effective spectral index can 
be inverted to obtain an estimate of the expected amount of curvature as a linear function 

of the spectral index F^: 

Tr — 2 

"'^^'^ ^ l0g{px/Pr) ^^^ 

where F^ is the effective spectral index in the radio band and pr is the fiux-weighted mean 
momentum of the electrons that produce the radio synchrotron emission. Although the 
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value oipx/Pr is unknown, px/Pr = [yxlvr^ , where vx and z/,. are the critical frequencies 
associated with px and Pr? respectively. Since the X-ray data used here span the range 
from 2 to 7 keV, vx ~ 9.0 x 10^^ Hz, the logarithmic mid-point of this band. Similarly, 
v^ ^ 660 MHz, the logarithmic mid-point of the range from 86 MHz to 5 GHz (Tabled]). 
These estimates yield log [pxIPr) = 4.57. This value is rather insensitive to the uncertainties 
in ux and z/^. For example, a change of a factor of 2 in the ratio fx/i^r corresponds to a 
change of only 3.3% in \og{px/Pr)- As shown in Figure [6] {dotted line), a relation of the 
form Oexp = (r^ — 2) /4.57 is consistent with the confidence contours for region 6. This 
result is remarkable because the contours in Figure [6] are sensitive to the relative X-ray and 
radio fluxes, while equation ([7]) has no such dependence. Had the relative normalizations 
been significantly different, it is unlikely that the contours in Figure [6] would have been 
consistent with the dotted line. If F^ is accurately represented by the best-fit index (i.e., if 
Pr ~ 1 GeV c~^) and if the cospatial radio fluxes are used, then a^xp = (2.221 — 2)/4.57 = 
0.048. This expected value is consistent with the fitted value (0.054 ± 0.006) at the 90% 
confidence level. If the peak radio fluxes are used, then Oexp = (2.198 — 2)/4.57 = 0.043, 
which is less consistent with the fitted value (0.033]'^q;oq8)- Like the expected radio flux, the 
expected amount of curvature lies between the cospatial and peak values. If the uncertainty 
associated with the radio flux leads to an uncertainty in the amount of curvature comparable 
to the statistical uncertainty, then the amount of curvature is most likely 0.05 ± 0.01. 

Not only is the best-fit amount of curvature consistent with equation ([7]), it is also 
consistent with the amount of curvature predicted for the proton (not electron) spectrum 
of SN 1006. For example. Figure [7| shows the best-fit electron spectra for region 6. The 
spectra are not curved below p = 1 GeV c~^, which was true during the fitting process, 
because the available data poorly constrain the shape of the electron spectrum in this mo- 
me ntum range. The top pair of dotted and dot-dashed curves are predictions described 



by lEUison et al.l (120001 ) for the proton (not electron) spectrum of SN 1006. These spectra 



have been normalized to match the solid black curve at E = 0.9 GeV. Since Ellison et al. 



(I2OOOI ) computed their spectrum assuming a spectral index F = 2.0, it is not appropriate to 
compare the solid curve with the upper pair of dotted and dot-dashed curves. However, the 
lower pair of dotted and dot-dashed curves have been multiplied by the energy- dependent 
factor [E/{0.9 GeV)]^°'^ to match the best-fit spectral index F = 2.2. Aside from a possible 
difference in the assumed value of the cutoff momentum, the lower pair of theoretical curves 
match the solid black curve remarkably well. 

The evidence of curvature is based on a set of assumptions. One assumption is that the 
electron spectrum has the form of equation (121). This functional form excludes the effects 



of a synchrotron cooling break (jVolk et al.ll2005l ). If the magnetic field strength is as large 



as 150 fiG, then synchrotron cooling may significantly steepen the highest-energy end of 
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the electron spectrum (JKsenofontov et al.ll2005l ). Since the fits are sensitive only to the net 
amount of curvature between radio and X-ray frequencies, inclusion of a cooling break would 
require a larger amount of curvature (i.e., more flattening than the amount of curvature 
reported here) to compensate for the spectral steepening of the cooling break. Another 
assumption is that the shape of the radio spectrum in each of the 13 regions is the same as 
the shape of the radio spectrum of the entire remnant. If this assumption is invalid, then the 
best-fit values of the spectral index and curvature may be inaccurate for any given region. 
Yet, the mean values of the index and curvature are insensitive to such spatial variations. 
A third assumption is that the shape of the synchrotron spectrum is uniform inside each 
of the 1 3 regions. Sorae you ng remnants, such as Cas A (jAnderson fc Rudnickl Il996l ) and 
Kepler (JDeLaney et al.ll2002l ) , exhibit evidence of radio spectral variations from one region to 
another. Therefore, the combination of the synchrotron spectra of several small-scale features 
might naturally yield a curved composite spectrum even if the synchrotron spectrum for each 
feature is just a power law. Small extraction regions were used to minimize this problem. 
A fourth assumption is that our estimates of the fraction of the 843 MHz flux from each 
region are accurate. As discussed above, the size of the beamwidth leads to some uncertainty 
about the proper radio fluxes to use for each region. Yet, even the worst-case scenario, the 
results obtained using the peak radio fluxes, requires some spectral curvature. High-quality 
radio and X-ray spectra for many, small, spatially resolved features in the remnant would 
eliminate the need for many of these assumptions. Unfortunately, these kinds of data are 
not available for SN 1006. 

Evidence of c urvature in the synchrotro n spectra of supernova remnants has been re- 
ported previously. [Reynolds fc EllisonI (119921 ) analyzed the remnant-integrated radio spectra 
of Kepler, Tycho, and SN 1006. They find hints of curvature for Tycho and Kepl er. Un- 
fortunately, a careful study of the radio emission from Kepler (JDeLaney et al.ll2002l ) reveals 
intrinsic variations in the radio spectral index from region to region. Therefore, it is not 
possible to dismiss the idea that the evidence of curvature in the remnant-integrated spec- 
trum is the result of the c ombination of several power-law spectra with different spectral 
indices. iJones et al.l (120031 ) minimized this problem for Cas A by analyzing the radio-to- 
infrared synchrotron spectra of small, selected features in the remnant. The results of their 
analysis indicate that the spectra of these regions flatten with increasing energy. We find 
that their spectra can be fitted with our curved spectral model if the curvature parameter 
a = 0.06 ± 0.01. It is interesting that this value is consistent with the mean amount of 
curvature for SN 1006 (0.05 ± 0.01). Based on an analy sis of X-ray data for a small feature 
on the northeastern shell of RCW 86, IVink et al.l (120061 ) report evidence of curvature in the 
synchrotron spectrum of this remnant. This claim would be strengthened if the data were 
fitted with a model that includes the spectral index as a free parameter and if there were 
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similar evidence for more than one feature in the remnant. 

The evidence of curvature suggests that the pressure exerted by nonthermal particles is 
large enough to modify the structure of the shock. The cosmic-ray pressure is most likely 
dominated by nuclei, not electrons. For example, analyses of the synchrotron spectrum of 
SN 1006 have been used to estimate the total energy of the nonthermal electrons in the 
remnant. These estimates, which depend on the assumed value of the mean magnetic field 
strength, include 5 x 10^^ ergs (for B = 150 fxG; Berezhko et al. 2002), 9 x 10"^^ ergs (for 
B = 40 /iG; Allen et al. 2001), and 7 x 10^^ ergs (for S = 3 /iG; Dyer e t al. 2001). The 



remnant was recently observed with the H.E.S.S. gamma-ray tele scope (lAharonian et al. 
20051 ) . If the upper limits on the TeV gamma-ray flux reported by lAharonian et al.l ( 20051 ) 
are accurate (cf. Tanimori et al. 1998), then a combination of these constraints and the 
measured synchrotr on flux can be used to s how that the mean magnetic field strength must 
be at least 25 ^G ( lAharonian et al.ll2005l ). In this case, the total energy in cosmic-ray 
electrons is probably no larger than about 10^^ ergs. Since this energy is a small fraction of 
the total internal energy (~ 10^^ ergs), the nonthermal electron pressure is most likely a small 
fraction of the ram pressure at the shock. Yet, the total energy in cosmic-ray nuclei, which 
may be about 2 orders of magnitude larger than the total energy in cosmic-ray electrons, 
could be large enough to modify the structure of the shock. 



3.2. Cutoff Frequency and Diffusion Coefficient 

In addition to providing evidence of curvature, fits to the synchrotron spectral data 
provide a measure of the frequency at which the synchrotron spectrum is cut off (eq. [S]). 
The results for the cutoff critical frequency are listed in Tables |2Hni A comparison of Ta- 
bles [2] and [3] (and of Tables H] and [5]) reveals that the frequencies for the curved spectra are 
somewhat lower than the frequencies for the uncurved spectra because positive curvature 
causes a synchrotron spectrum to flatten with increasing energy. To compensate, the X-ray 
synchrotron spectrum is steepened by reducing the cutoff frequency. The solid and dashed 
black lines in Figure [7] show that the two effects more or less offset one another to produce 
electron (and, hence, synchrotron) spectra that have nearly the same shape in the X-ray 
synchrotron-emitting band. 



Bamba et al.l (l2003r ). iRothenfiug et al.l (J200J), and iBamba et al.l (120081 ) fitted X-ray 



spectra for the bright rims of SN 1006 with models that include the XSPEC synchrotron 
component srcut. Since srcut does not include spectral curvature, their results should 
be compared with the r esults listed in Table 3. To the extent that the regions used by 
Bamba et al.l (|2003l ) and IRothenfiug et al.l (|2004l ) are similar to the regions used here, this 
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comparison reveals that their cutoff frequencies are systematically higher than the frequencies 
listed in Table [31 This difference is due, in part, to a difference in the spectral indices. Our 
best-fit spectral indices (F = 1.98-2.12; Table [3]) are lower than the indices used by Bamba 
et al. (2003; T = 2.14) and Rothenflug et al. (2004; T = 2.2). While an index of about 2.2 
is appropriate for the radio synchrotron-emitting electrons (i.e., the electrons at momenta 
p ~ 1 GeV c~^), an index of 2.0 seems to better describe the X-ray synchrotron-emitting 
electrons (i.e., the electrons at momenta p ~ 10 TeV c~^). For example, if the electron 
spectrum is not curved, then the mean index is 2.031 ± 0.017 (Table [3]). If curvature is 
included, then the effective mean index at 10 TeV is 2.005 ± 0.027 [i.e., 2.221 ± 0.013 - 
(0.054 ±0.006) log(10 TeV/1 GeV); Table [2]. In both cases, the results imply that the index 
near the cutoff of the electron spectrum is about 2.0, which is consistent with the model of 



Berezhko et al.l (|2002| ). Berezhko et al. note that the electron spectrum is steeper at radio 
synchrotron-emitting momenta than at X-ray synchrotron-emitting momenta because GeV 
electrons sample only a portion of the shock while the highest-energy electrons "see" the full 



strength of the shock. Although the same argument applies to the results of iBamba et al. 



(120081 ). who use a spectral index F = 2.14, their best-fit break frequency for the entire 
northeastern rim is lower than the mean value of i>m in TableO [i^m = (9.07±0.98) x 10^^ Hz]. 
If the normalization of the synchrotron spectrum is too high, then the best-fit frequency may 
be reduced. While the normalization reported by Bamba et al. (2008; 7.72 Jy at 1 GHz) is 
approximately correct for the entire northeastern rim, only the flux density of the portion of 
the rim associated with the very narrow region in which the X-ray synchrotron emission is 
produced should be used. For these reasons, we beheve that the cut-off frequencies reported 
here are relatively accurate. 

The frequencies listed in Tables [2] and [3] are plotted in Figure [8] as a function of the 
position angles of the regions. This angle is measured relative to the location of the center 
of the remnant (Fig. [T]) and is expressed in degrees from north (0°) through east (90°). As 
shown in Figure [H], there may be an azimuthal variation in the value of the cutoff frequency. 
A constant value [i.e., the weighted mean value of (4.98 ± 0.67) x 10^^ Hz, 90% con f idence 



level uncertainties] can be excluded at the 2.5 a confidence level. iRothenflug et al.l (120041 ) 
also report evidence of an azimuthal variation in the cutoff frequency. Since the frequency 
Um oc p'^B, these results imply that the cutoff momentum of the electron spectrum and/or 
the magnetic field strength varies along the northeastern rim. Since it is not possible to 
independently determine both pm and B using the synchrotron spectral data alone, the 
present analysis does not provide a means of identifying the cause of a variation in the cutoff 
frequency. 

As described in the Appendix, it is possible to use a fitted value of the cutoff critical 
frequency to constrain the mean electron diffusion coefficient k. The upper limit on k, as a 
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fraction of the mean Bohm diffusion coefficient kb, depends on the forward-shock velocity 
Ui, the cutoff critical frequency Um, and a function / of the compression ratio r and the 
ratio of the upstream t o downstream magnetic field strengths B^/B-} (eas. | A5|- | A8|N^| If 
u-\ = 2300-5000 km s~^ (JLaming et al.lll996l : iDwarkadas fc Chevalierlll998l : iGhavamian et al. 
2002h . Um = 1.1 X 10^^ Hz (Table [2]), and / = 0.1875 (i.e., r = 4 and S1/S2 « 0), then 



K < (4.5-21)kb. If the compression ratio is 5.2 (JKsenofontov et al.ll2005l ) instead of 4, then 
f = 0.155 and R < (3.7-18)fi;B. If the cutoff frequency is about 10^® Hz (JRothenflug et al. 



20041 ) instead of 10^^ Hz, then R < (0.5-2.3)/t;b. In each of these cases, the highest-energy 
electrons have diffusion coefficients nearly as small as the Bohm coefficient (i.e., are being 
accelerated about as fast as possible). 



Similar results for the diffus ion coefficient are reported by lBamba et al.l(l2003l ). lYamazaki et al. 



( 120041 ). and lParizot et al.l (120061 ) . The latter two of these three use equations similar to equa- 
tion (]A5P but obtain significantly smaller limits for the diffusion coefficient because they use 
much larger cu toff critical frequenc ies (9.1 x 10^^ and 2.5 x 10^^ Hz, respectively, instead of 
1.1 X lO'^^ Hz). iBamba et al.l (120031 ) take a different approach. They assume that the width 
{ipid) of the part of a filament that is in the upstream region is determined by the diffusion 
length [ki/ui) in this region. In this case. 
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(9) 



where ki, k^i [= pc/{3eBi)], Bi, I'm,!, and ipi are the diffusion coefficient, Bohm diffu- 
sion coefficient, mean magnetic field strength, cutoff critical frequency, and filament width jj 
respectively, in the upstream region and Mi is the velocity of the forward shock. The nu- 
merical constant in equation (j9]) was computed assuming that the value of sin^i is 71/4, 
the isotropic mean value. Like equation (]A5p . equation (Q is valid only at the momen- 
tum p = pm- If the value of Bi is comparable to a typical interstellar field strength (e.g., 
3 /iG), then, for reasonable values oiipi, d, Ui, and i^m,!, 1^1 is uncomfortably small compared 
with the Bohm diffusion coefficient. There are at least three explanations for this appar- 
ent dilemma. One explanation is that, near the shock, Bi may be amplified by cosmic-ray 



'"^The corresponding upper limit on the diffusion length Im w k/ui (JLagage fc Cesarskvlll983l ). Using 
eqs. O and (1X41) yields U oc fmv' H^Bo^^'^ . If / = 0.1875, ui = 2300-5000 km s'^, Vm,2 = lO" Hz, and 
B2 > 25 fiG JAharonian et al.fboosl ). then /,„ < 0.17-0.36 pc (i.e., U/d < 16"-34" if d = 2.2 kpc; Winkler 
et al. 2003). 



^Here ipi is the actual angular width, not the observed angular width. 
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streaming flLucek fc Bel]|l200d : Bell fc Lucekll200ll : lBel]|l2004l ). For example, equation (15) of 
Bell fc Lucekl (120011 ) yields Bi = 30 /iG (see Ksenofontov et al. 2005) and, hence, ki ~ /€b,i, 
if the upstream mass density pi = 2.4 x 10~^^ g cm~^ (i.e., rii = 0.1 cm~^), the shock 
velocity Ui = 2300-5000 km s~^, and the cosmic-ray pressure is abo ut 30%-14% of the ram 
pressure, respectively. Another explanation, which is mentioned by lYamazaki et al.l (120041 ). 
is that the rate of acceleration can exceed the Bohm limit (i.e., ki can be less than kb,i) if 
the magnetic field is parallel to the shock (see Jokipii 1987). A third explanation is that the 
width of a filame nt represents the size of the region in which the magnetic field is strong 
( iPohl et al.ll2005l ) instead of the diffusion length. In this case, equations ([8]) and (Q are not 
useful because an assumption upon which they are based is not valid. 

Independent of the technique used to estimate or constrain the diffusion length, it seems 
likely that the highest-energy electrons in SN 1006 are diffusing close to the Bohm limit. 
Otherwise, the rate of acceleration is too low for electrons to reach energies high enough 
to produce the observed X-ray synchrotron radiation. Yet, it is important to note that the 
limit on the diffusion coefficient described in the Appendix (and the value obtained using 
eqs. [8] and [9]) is only an estimate, for the following reasons: 

1. The fitted values of the cutoff frequency are sensitive to the assumed shape of the 
electron spectrum at momenta near the cutoff in the electron spectrum. The actual functional 
form of the electron spectrum may be significantly more complicated than equation ([2]). As 
an example of the extent to which differences in the shape of the electron spectrum can affect 
the limit, consider the results obtained using the models with and without curvature. If the 
curved model is used, then the maximum cut-off critical frequency is 1.1 x 10^^ Hz (Table [2]) 
and R < (4.5-21)kb (for / = 0.1875, B1/B2 ^ 0, and ui = 2300-5000 km s'^). If the model 
does not include curvature, then z/^ = 3.0 x 10^^ Hz (Table [3]) and R < (1.7-7.8)kb. 

2. There is evidence that the cutoff frequency varies from one region to another in 
SN 1006 (Fig. [HI Rothenfiug et al. 2004). Here the frequency used to compute the limit is 
the largest value in Table [2l Therefore, our constraint on the diffusion coefficient of SN 1006 
applies only to the region where the cutoff frequency is largest. The limits for other regions 
of the remnant are larger. 

3. The limit is computed assuming that the electron momentum p = pm- Therefore, the 
limit only applies to electrons at the cutoff momentum [i.e., R, = R{p = Pm)]- An advantage 
to using this momentum is that the limit does not depend on the functional form of R{p). 

4. The limit is based on the assumption that the cutoff is due to synchrotron losses. 
If some other process, such as the escape of electrons from the acceleration region or the 
length of time over which particles are accelerated, determines the cutoff momentum, then 
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the diffusion coefficient is lower than the right-hand side of equation flA5|l . If this equation 
yields a limit close to 1, then the results indicate that synchrotron losses are important for 
electrons that have momenta p > Pm- 

5. The value of the function / depends on the unknown compression ratio and the 
unknown ratio of the upstream to downstream magnetic field strengths (eq. [A7] ). A value 
of / = 0.1875 (i.e., r = 4 an d B^ /By, rj 0) is used here because the compression ratio is prob 
ably at least as large as 4 (IBlondin &: EUisoru 1200 ll : iKsenofontov et al.l l2005l : IWarren et al. 



20051 : ICassam-Chenai et al.ll2008l ). If r > 4, then / < 0.1875 and our upper limit on the 



diffusion coefficient is overestimated. Likewise, the upper limit is overestimated if B1/B2 is 
significantly larger than zero. 

6. Equation flA5p was derived assuming an isotropic pitch-angle distribution. While 
such a distribution is expected if electrons diffuse in the Bohm limit, the actual distribution 
may be different. 



Conclusions 



We have performed a joint spectral analysis of some Chandra ACIS X-ray data and 
MOST radio data for 13 small regions along the bright northeastern rim of the supernova 
remnant SN 1006. The data were fitted with a model that includes a synchrotron emission 
component. This component is based on an electron spectrum that has a momentum- 
dependent spectral index. The rate of change in the index for each decade in momentum is 
a free parameter of the fit. If the assumptions described in § 13.11 are valid, then the results 
of the spectral analysis, which show that the synchrotron spectra of SN 1006 are curved, 
can be interpreted as evidence of curvature in the GeV-to-TeV electron spectra. The mean 
amount of curvature in the electron spectra is qualitatively consistent with pr edictions of 
the amount of curvature in the proton spectrum of SN 1006 (JEUison et al.ll2000l ). The best- 



fit power-law index at 1 GeV (i.e., at radio synchrotron-emitting momenta) is 2.221;!^q'qJ2- 
Including the effect of curvature, the effective spectral index at about 10 TeV (i.e., at X-ray 
synchrotron-emitting momenta) is 2.005 ± 0.027 ( 90% confidence level u ncertainties). This 
effective index is consistent with the predictions of iBerezhko et al.l (120021 ). 



The evidence of curved electron spectra suggests that cosmic rays are not "test" par- 
ticles. The cosmic-ray pressure at the shock is large enough to modify the structure of the 
shock. Since nonthermal electrons contain only about 0.1% (i.e., 10^^ ergs) or less of the 
total internal energy, the results provide indirect evidence of a much more energetic popu- 
lation of cosmic-ray protons. Collectively, the evidence of (1) spectral curvature in Cas A 
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fl Jones et al.ll2003[ ). RCW 86 flVink et al.ll2006h and SN 1006, (2) an unusually low electron 
temperature in I E 0102.2— 7219 (IHughes et al.l 120001 ), and (3) a compression rat io greater 
than 4 in Tycho (j Warren et al.l |2005| ) and SN 1006 (jCassam-Chenai et al.l |2008| ) suggests 
that efficient particle acceleration may be a common feature of young, shell-type supernova 
remnants. 



The results of the spectral analysis also determine the "cutoff critical frequency" Vm- 
This frequency seems to vary from one region to another (see also Rothenfiug et al. 2004), 
which implies that the exponential cutoff momentum of the electron spectrum and/or the 
strength of the magnetic field varies. It is not possible to identify the cause of the variation 
using the synchrotron spectral data alone. 

As described in the Appendix, the cutoff frequency can be used to set an upper limit 
on the mean diffusion coefficient R of the highest-energy electrons. Aside from the cutoff 
frequency, this limit depends (strongly) on the velocity of the forward shock ui and (weakly) 
on the compression ratio r and the ratio of the upstream to downstream magnetic field 



strengths B1/B2. If i^m = 1-1 x 10^^ Hz, Ui = 2300-5000 km s" 



4, and B1/B2 ^ 0, 



then K < (4.5-21)/t;b, where Rb is the mean Bohm diffusion coefficient. This result implies 
that at least some of the highest-energy electrons in SN 1006 diffuse close to the Bohm limit 
(i.e., are accelerated about as fast as possible), which provides additional support for the 
idea that Galactic cosmic rays are predominantly accelerated by the shocks of supernova 
remnants. 
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A. Diffusion coefRcient 



This appendix describes how measurements of or inferences about the shock velocity 
and cutoff frequency can be used to place an upper limit on the electron diffusion coefficient. 
Since the mean rate of synchrotron losses cannot exceed the mean rate of energy gains at 
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momenta below the cutoff of the electron spectrum, 



\ / ace \ / sync 

in this range. The rate of acceleration at p = p^ ^ mc is given by 

(Lagage & Cesarsky 1983). If the rates of synchrotron losses in the upstream and downstream 
regions are weighted by the mean residence times in these regions [Ati = Aki/^uiv) and 
At2 = Ak2/{u2v)] Webb et al. 1984], then 

dE\ e^ 2 ('^It32 ■ 2q , f^2r^2 ■ 2n\ ff^l , '«2V^ f \ Q\ 

Pm — ^1 sm 6*1 H B^ sm^ 6*2 \ (A3) 



'^^/sync 67160/71"^ C^ ™ VMl U2 J \Ul U2 

(Blumenthal & Gould 1970), where ki and k,2 are the upstream and downstream electron 
diffusion coefficients perpendicular to the shock, ui and U2 are the speeds of the upstream 
and downstream material relative to the shock, v is the velocity of a particle, eo is the 
permittivity of free space, Bi and B2 are the mean upstream and downstream magnetic field 
strengths, and 9i and 62 are the pitch angles between the electron momentum and magnetic 
field vectors in the upstream and downstream regions, respectively. If sin 9 in equation (jlj) 
is replaced by the isotropic mean value of 7r/4 and if sin^ 61 and sin^ 62 in equation (1A3P are 
replaced by the isotropic mean value of |, then a combination of equations (jlj), (jAip . (]A2p . 
and dSH]) yields 



K 27Tieomc fu\ 

Kb ~ 16e2 i/„o 



< l!i!i^A^ (A4) 



-1 



where 



^ ''■^^*'la45)(liF^)^(u;^) • (^^) 



' r2 r2 \ / r2 r2~ 

..= (^^, + ^.2^^ + ^) , (A6) 

R^ = ('^/.B,! + ^KB,2l (^ + ^] \ (AT) 



Ml U2 / \ Ui U2 



Ml M2 / V^l U2 

the upstream and downstream Bohm diffusion coefficients kb.i and fi;B,2 are equal topmc/(3e-Bi) 
and pm.c/(3ei?2), respectively, at a momentum p = pm, i^m,2 is the critical frequency of an 
electron with this momentum in a magnetic field B = B2, 

r — 1 

^ = r[r + {B,lB2)Y ^^^'^ 
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and the compression ratio r = uxju^. The value of 0.1875 in equation flASP was calculated 
assuming r = 4 and BxjB^ ~ 0. If the compression ratio is greater than 4 or B^jB^ is 
significantly larger than zero, then / < 0.1875. In no case can / be l arger than .25. I f 
Bx = B2, then equation flA5p is identical to t he analogous equatioii i n IStage et al.l (120061 ). 
E quation (IA5D is simila r to equation (22) of lAharonian fc AtovanI (119991 ) . equation (12) 
of iLazendic et al.l (|2004[) . equation (A. 4) of lYamazaki et al.l (120041 ). and equation (22) of 



Parizot et al. 



(l2006l ). except these authors use the peak frequency instead of the critical 



frequency. The latter three also use sin 6* = 1 instead of 7r/4. 
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Fig. 1. — Color-coded ACTS image of the northeastern rim of SN 1006. Red, green, and blue 
correspond to the energy bands 0.4-1, 1-2, and 2-7 keV, respectively. Most of the emission 
at energies greater than 1 keV is concentrated in the whitish filaments along the rim. The 
edges of the six square CCDs used to observe the remnant are evident. The relatively high 
brightnesses of the CCD with the circle on it and the CCD in the lower left corner are due to 
enhanced low-energy sensitivities for these two devices. The circle at a = 15^03™51.56** and 
S = — 41°51'18.8" (J2000) indicates the location of the nominal aim point of the telescope. 
The asterisk at a = 15^02™51.7*^ and S = — 41°56'33.0" is the location of the center of the 



supernova remnant fj Winkler &: Longlll997l ). The 13 squares are the X-ray spectral extraction 
regions. 
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Fig. 2. — Radio spectrum for the entire supernova remnant SN 1006. Top: The data 
points (Table [1]) and three models. The dashed line is the best-fit power law, 5'(z/) = 
17.6^4g[z//(l GHz)]~" Jy, where a = O.GOlg'Qg. The dotted line is a power law using Green's 
(2006) spectral parameters [S^u) = 19[z//(l GHz)]^°'® Jy]. The solid line is the best-fit 
(curved) synchrotron model for region 6 divided by the factor 0.00175 (see Table [2]). Bottom: 
Differences between the data points and the solid line, divided by the uncertainties in the 
data points. 
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Fig. 3.— An 843 MHz MOST image of SN 1006, courtesy of R. Roger jRoeer et alJIlOSsh . 
The ellipse at lower right indicates the half-power beamwidth. The 13 black squares are the 
locations of the X-ray spectral extraction regions. These regions were used to obtain the 
"cospatial" radio fluxes (Tables [2] and [3]). The red squares are the locations used to obtain 
the peak radio fluxes (Tables H] and [5]). The ratios of the flux densities of the black and 
red regions to the total flux density are listed in the ( columns of Tables [2H51 The bright 
spot at about a = 15^04™ 04'^ and 6 = — 41°55'48" is produced by an extragalactic source 



(JReynolds fc Gilmorelll986l ). 
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Fig. 4. — ACIS spectrum for region 6 (see Fig. [T]) of the bright northeastern rim of SN 1006. 
Top: Sum of the source and background spectra (data points) and the sum of the best- 
fit model and background spectra {solid line). The model displayed here includes spectral 
curvature as a free parameter (see Table [2]). Bottom: Differences between the data points 
and the solid line, divided by the uncertainties in the data points. 
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Fig. 5. — X-ray emission profiles along a 49"- wide strip passing through region 6 toward 
the center of SN 1006. The black histogram in the top panel is the 2-7 keV ACIS data. For 
comparison, the black line in the bottom panel is a model based on the assumption that the 
emission is produced in a uniformly emitting shell that is 6" wide. The red curves are the 
data and model smoothed to the resolution of the MOST image. (The horizontal red bars 
are the 1 a widths of the Gaussian smoothing function.) The reduction in spatial resolution 
causes the peak to shift downstream, but the offset is small compared with the width of 
the extraction region. (The blue lines are the boundaries of region 6.) Therefore, it seems 
unlikely that the offsets between the X-ray and radio peaks (see Fig. [3] and Tables H] and [5]) 
are due to the different X-ray and radio resolutions. 
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Fig. 6. — The 1, 2, and 3 a confidence contours for region 6 (see Fig. [T]) in tlie parameter 
space defined by the electron spectral index F and electron spectral curvature a (see eq. [2]). 
The sohd black and dashed red contours are the results obtained using the cospatial (Table [2]) 
and peak (Table Hj) radio fluxes, respectively. The plus signs indicate the best-fit values of 
the index and curvature. The dotted black line is the expected relationship between F and 
a (eq. [7]). The dashed black line is the line along which the electron spectrum is not curved 
(i.e., a = 0). For region 6, a positive curvature (i.e., a result above the dashed line) is favored 
at about the 2.7 and 1.6 a confidence levels for the black and red contours, respectively. 
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Fig. 7. — Best-fit electron number density spectra for region 6 (see Fig. [T]). The solid and 
dashed black lines are the curved (Table [2]) and uncurved (Table |3]) models, respectively, 
using the cospatial radio flux. The solid red line is the curved model using the peak radio 
flux (Table H]). The upper dotted and dot-das hed curves ar e the m odels for the proton (not 
electron) spectrum of SN 1006 presented by lEUison et al.l (120001 ). These lines are plotted 
only for momenta p > mc (i.e., the range of momenta to which our flts are sensitive) and 
are normalized to the solid black line at a kinetic energy of 0.9 GeV. As described in the 
text, the lower dotted and dot-dashed curves are the same pair of models multiplied by 
[E/{0.9 GeV)]~°'^. The amounts of spectral curvature in the solid lines are consistent with 
the amounts of curvature in the lower dotted and dot-dashed curves. From left to right, the 
three bracketed energy bands contain the electrons that are primarily responsible for the 
observed thermal bremsstrahlung, radio synchrotron, and X-ray synchrotron emission. 
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Fig. 8. — Cutoff frequency as a function of position angle (measured counterclockwise 
from north through east). From right to left, the 13 pairs of points are for regions 1-13, 
respectively (see Fig. [1]). The black data points and 90% confidence intervals are the results 
obtained using a curved electron spectrum and the cospatial radio fluxes (Table [2]). The red 
data points, which have similar confidence intervals, are the results obtained using the peak 
radio fluxes (Table HI). The dashed line is the weighted mean value (4.98 x 10^^ Hz) of the 
black points. The corresponding 90% confidence level interval [(4.31-5.65) x 10^^ Hz] lies 
between the two dotted lines. A constant value for the cutoff frequency can be excluded at 
the 2.5 (J confidence level. 
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Table 1. Radio Data 



V 


S 


A5 




(GHz) 


(Jy) 


(Jy) 


Refs. 


0.086 


77 


7.7^ 


1 


0.16 


46 


13.8 


1 


0.408 


32.3 


3.2 


1 


0.408 


30.8 


3.P 


2 


0.635 


25.5 


3.8 


1 


0.843 


17.5 


1.5 


3 


1.41 


15.2 


2.3 


1 


2.7 


10.4 


1.6 


1 


2.7 


9.9 


1.0^ 


1,4 


5.0 


i A 


1.2 


1 


5.0 


7.0 


0.7^ 


5 


5.0 


8.1 


0.8^ 


6 



Note. — The table lists the 
frequencies z/, flux densities 5", 
and 1 o flux density uncertain- 
ties AS* used for the spectral 
analysis. The quantities S and 
A^* are for the emission from 
the entire supernova remnant. 

'^Since no uncertainty was 
reported, it is assumed to be 
10% of S. The results of our 
spectral flts are insensitive to 
the assumed fractional uncer- 
tainty at least for fractions be- 
tween 5% and 20%. 

References. — (1) Milne 



-35- 

1971; (2) Stephenson et al. 
1977; (3) Roger et al. 1988; 

(4) Gardner & Milne 1965; 

(5) Kundu 1970; (6) Milne & 
Dickel 1975. 



Table 2. Best-Fit Parameters Using the Curved Model and "Cospatial" Radio Fluxes 





ax 


Sx 





c 






t'm 






Region 


(J2000) 


(J2000) 


(deg) 


(xlO-4) 


pa 


a 


(10^6 Hz) 


xVdof 


^Ax^' 


1 


15 04 10.7 


-42 00 15 


104.1 


10 


r, o .+0.16 

^•^^-0.16 


U.UOU_o.o59 


2.01?:^ 


8.4/18 


0.17 


2 


15 04 13.3 


-41 57 51 


94.8 


8.9 


9 97+0.13 


0.0711°:°^^ 


1-81?:^ 


12.9/22 


0.012 


3 


15 04 13.7 


-41 56 27 


89.6 


6.9 


9 97+0.13 
^•^'-0.13 


0.0801°:°^^ 


2.11- 


11.0/27 


0.00055 


4 


15 04 09.7 


-41 54 54 


83.5 


24 


9 99+0.12 


0.0511°:°^° 


1+2.1 
<J--'--1.2 


50.2/46 


0.0073 


5 


15 04 09.7 


-41 53 41 


78.8 


15 


9 90+0.11 
^•^'J-0.12 


0.060l°r7 


4.31- 


71.4/63 


0.00033 


6 


15 04 06.5 


-41 51 25 


69.7 


18 


9 91+0.11 
^•^-L-O.ll 


o.o5il°:°i 


5.81?:^ 


91.5/84 


0.0014 


7 


15 04 04.8 


-41 49 35 


62.8 


11 


2.191°- 


0.0361°-°^^ 


8.911^ 


69.2/60 


0.030 


8 


15 04 01.1 


-41 49 22 


60.8 


17 


2.201°:}! 


0.042l°-°26 


6.111^ 


69.4/71 


0.011 


9 


15 03 57.3 


-41 48 08 


55.4 


19 


2.201°:}} 


0.042l°:°f, 


6.01?:^ 


73.0/72 


0.010 


10 


15 03 52.6 


-41 46 58 


49.7 


17 


2.201°:}^ 


0.042l°-°26 


6.5+I2' 


71.4/69 


0.011 


11 


15 03 57.7 


-41 45 28 


47.9 


3.8 


2.201°:}} 


0.0571°:°! 


I0.8tl%' 


65.3/61 


0.00079 


12 


15 03 52.3 


-41 43 31 


40.8 


3.0 


9 91+0.11 
^•^J^-0.11 


o.o6oi°:°i 


8.2tl'f 


32.2/37 


0.0018 


13 


15 03 29.5 


-41 42 03 


25.8 


9.5 


9 9P:+0.11 
^•^3-0.12 


0.0701°:°^} 


q q+2.5 


35.5/42 


0.00023 


Mean 










9 991+0.013 
Z.ZZ1_QQ^2 


0.054 ±0.006 









00 



Notes. — For each region, the table includes a list of the right ascension (ax), declination (5x), and position- 
angle (0) coordinates of the center of the X-ray region, the 843 MHz flux divided by the total 843 MHz flux 
{(), the best-fit differential electron spectral index F, curvature parameter a, and cutoff frequency z/m, the value 
of x^ P^i' degree of freedom (x^/dof), and the probability (-PaxO ^^at the difference between the value of x^ 
in Table [3] and the value in this table is due to chance. The uncertainties, which include only the statistical 
contributions, are reported at the 90% confidence level. These uncertainties were used to compute the weighted 



mean values of F and a. Units of right ascension are hours, minutes, and seconds, and units of dechnation are 
degrees, arcminutes, and arcseconds. 

'^As described in the text, the spectral indices are not independent from region to region. 

^Including curvature as a free parameter significantly improves the quality of the fit for most regions. 
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Table 3. Best-Fit Parameters Using the Uncurved Model and "Cospatial" Radio Fluxes 





ax 


Sx 





C 






l^m 




Region 


(J2000) 


(J2000) 


(deg) 


(xlO-4) 


F 


a"" 


(10^6 Hz) 


xVdof 


1 


15 04 10.7 


-42 00 15 


104.1 


10 


9 1 9+0.08 
^•-"-^-0.08 





5.8+^i 


10.3/19 


2 


15 04 13.3 


-41 57 51 


94.8 


8.9 


9 08+0-07 





Q.strs 


19.2/23 


3 


15 04 13.7 


-41 56 27 


89.6 


6.9 


'^•"-''-'-0.07 





n v+6.1 
'•'-3.0 


22.9/28 


4 


15 04 09.7 


-41 54 54 


83.5 


24 


9 04+0-05 
-^•U4_o,o6 





c 0+3.1 


57.4/47 


5 


15 04 09.7 


-41 53 41 


78.8 


15 


1 08+0-04 





l&t\\ 


84.3/64 


6 


15 04 06.5 


-41 51 25 


69.7 


18 


1 qq+0.04 





9.21^-6 


101.7/85 


7 


15 04 04.8 


-41 49 35 


62.8 


11 


o nei+0.05 





14.7+^-3 


73.9/61 


8 


15 04 01.1 


-41 49 22 


60.8 


17 


9 nQ+0.04 
^•'-''-'-0.04 





9.8^to' 


75.9/72 


9 


15 03 57.3 


-41 48 08 


55.4 


19 


o no+0.04 
^•'-''-'-0.04 





9.4^^3 


79.6/73 


10 


15 03 52.6 


-41 46 58 


49.7 


17 


9 nQ+004 
^•'-''-'-0.04 





10.3+^:? 


77.9/70 


11 


15 03 57.7 


-41 45 28 


47.9 


3.8 


9 00+0-04 





28.4+?^-? 


76.5/62 


12 


15 03 52.3 


-41 43 31 


40.8 


3.0 


9 no+0.06 
^■'-'^-0.06 





29.6+?^i 


41.9/38 


13 


15 03 29.5 


-41 42 03 


25.8 


9.5 


9 nn+005 





8.811? 


49.1/43 


Mean 










9 n'^i +0-017 










Notes. — See the caption for Table O Units of right ascension are hours, minutes, and 
seconds, and units of declination are degrees, arcminutes, and arcseconds. 

'^The curvature parameter was fixed at zero. 



Table 4. Best-Fit Parameters Using the Curved Model and Peak Radio Fluxes 



Region 



an 6r A^ 

(J2000) (J2000) (arcmin) 



c 

x 10-^1 



(1016 Hz) xVdof 






1 


15 04 07.1 


-42 00 05 


0.68 


17 


9 9-1+0.16 
^•^-■--0.16 


0.035+°:°^^ 


2 1+^-^ 


8.5/18 


0.35 


2 


15 04 07.7 


-41 57 44 


1.04 


22 


9 99+0.14 
^•^^-0.15 


0.042+°:°^? 


2.0+f^ 


13.1/22 


0.15 


3 


15 04 04.4 


-41 56 25 


1.73 


30 


9 9r.+0.12 


U.UO±_o.o38 


2.6+?:° 


11.3/27 


0.18 


4 


15 04 04.6 


-41 55 00 


0.95 


30 


9 91+0.12 


0.043+°:°^J 


3.2+li 


50.2/46 


0.021 


5 


15 04 06.1 


-41 53 48 


0.67 


22 


9 91+0.11 
^•^-■--0.12 


OMSt'oTr 


A.Qtli 


71.5/63 


0.0037 


6 


15 04 02.6 


-41 51 40 


0.77 


31 


9 iq+0.11 

^■-'-y-0.12 


U.UOO_oo26 


6.61IS 


91.8/84 


0.040 


7 


15 03 59.9 


-41 50 02 


1.03 


26 


2 16+°-ii 

^•-LO-O.ll 


nnin+o.027 

U.UiU_o.o28 


10.9+^-^ 


69.6/60 


0.55 


8 


15 03 58.6 


-41 49 37 


0.53 


20 


2 20+°-ii 


Q Q07+O.O27 


6.3+f? 


69.5/71 


0.023 


9 


15 03 54.1 


-41 48 31 


0.72 


22 


9 9n+o-ii 

^•^'-'-0.12 


0.037+°:°^? 


6.2+1^ 


73.0/72 


0.023 


10 


15 03 50.5 


-41 47 18 


0.52 


19 


2 20+°-ii 
^■^'-'-0.12 


0.038l°:°i 


6.611^ 


71.4/69 


0.019 


11 


15 03 47.8 


-41 46 56 


2.35 


16 


2.16:^°:^} 


o.oi3+°:°i 


16.0+^^6^ 


65.7/61 


0.46 


12 


15 03 41.7 


-41 45 48 


3.02 


15 


2.iQtrA 


OMOt'o'ofi 


i2.o:^^y 


32.6/37 


0.59 


13 


15 03 25.2 


-41 43 45 


1.88 


17 


9 99+0.12 
^•^^-0.12 


U.UO/_o.031 


q 7+2.8 
<J-'-1.5 


35.6/42 


0.0062 


Mean 










2.198 ±0.011 


U.UOO_o.o08 









00 
CO 



Notes. — For each region, the table includes a list of the right ascension (aji) and declination {6fi) coordinates 
of the center of the radio region, the angular distance between the centers of the X-ray (see Table [2]) and radio 
regions A\E', the 843 MHz flux divided by the total 843 MHz flux {(), the best-fit differential electron spectral 
index F, curvature parameter a, and cutoff frequency z/^, the value of x^ P^r degree of freedom (x^/dof), and 
the probability (-PaxO fhat the difference between the value of x^ i^i Tableland the value in this table is due 
to chance. The uncertainties, which include only the statistical contributions, are reported at the 90% confidence 



level. These uncertainties were used to compute the weighted mean values of F and a. Units of right ascension are 
hours, minutes, and seconds, and units of declination are degrees, arcminutes, and arcseconds. 

^As described in the text, the spectral indices are not independent from region to region. 

'^Including curvature as a free parameter improves the quality of the fit for every region. 
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Table 5. Best-Fit Parameters Using the Uncurved Model and Peak Radio Fluxes 





aR 


5r 


A^ 


C 






T^m 




Region 


(J2000) 


(J2000) 


(arcmin) 


(xlO-4) 


r 


a"" 


{W Hz) 


xVdof 


1 


15 04 07.1 


-42 00 05 


0.68 


17 


9 1 q+0.08 

^•-'^'-'-0.08 





4.3+f? 


9.4/19 


2 


15 04 07.7 


-41 57 44 


1.04 


22 


9 1 1 +0.07 

^•-"^-"--o.os 





4.3:^?i 


15.2/23 


3 


15 04 04.4 


-41 56 25 


1.73 


30 


9 1 Q+0.07 





A2tll 


13.1/28 


4 


15 04 04.6 


-41 55 00 


0.95 


30 


o ne;+0-06 
^•"J^-0.05 





5 5+2-9 


55.6/47 


5 


15 04 06.1 


-41 53 48 


0.67 


22 


2 01+005 





7.3^|S 


79.9/64 


6 


15 04 02.6 


-41 51 40 


0.77 


31 


o r,e-+0.04 





8.9t|^ 


96.1/85 


7 


15 03 59.9 


-41 50 02 


1.03 


26 


9 19+0.05 
^•-^^-0.04 





i2.6:!:I:I 


70.0/61 


8 


15 03 58.6 


-41 49 37 


0.53 


20 


2Mtroi 





9.611^ 


74.6/72 


9 


15 03 54.1 


-41 48 31 


0.72 


22 


2.04tHl 





9.211^ 


78.2/73 


10 


15 03 50.5 


-41 47 18 


0.52 


19 


9 n4+0.04 
^■'J^-0.04 





io.i:^IS 


76.9/70 


11 


15 03 47.8 


-41 46 56 


2.35 


16 


9 1 9+0.04 
^•-"-^-0.05 





2o.o:^^V 


66.2/62 


12 


15 03 41.7 


-41 45 48 


3.02 


15 


9 1 q+0.06 

^•-'-'-'-0.06 





15.0;^^4^ 


32.9/38 


13 


15 03 25.2 


-41 43 45 


1.88 


17 


9 r,4+0.05 





7.511^ 


43.1/43 


Mean 










9 n7Q+0.021 
Z.U(0_o.020 










Note. — See the caption for TableHl Units of right ascension are hours, minutes, and seconds, 
and units of declination are degrees, arcminutes, and arcseconds. 

^The curvature parameter was fixed at zero. 



